The right tail exponent of the Tracy- Widom-/3 

distribution 



Laure Dumaz Balint Virag 

January 13, 2013 



Abstract 

The Tracy- Widom j3 distribution is the large dimensional limit of the top eigenvalue 
of (3 random matrix ensembles. We use the stochastic Airy operator representation to 
show that as a — )• oo the tail of the Tracy Widom distribution satisfies 

P {TWp >a) = a-i/^+°(i) exp (-\^cc'l^ 





Flowlines for the ODE of X - -^B 



1 Introduction 



For /3 > fixed, we examine the probability density of Ai > A2 > • • • > A„ G M given by : 



in which Z„ is a normalizing constant. This family of distribution is called the /3-ensemble. 
When P = 1,2 or 4, this is the joint density of eigenvalues for respectively the Gaussian 
orthogonal, unitary, or symplectic ensembles of random matrix theory. But the law ([T]) has a 
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ph ysical sense for the as i t descr ibes a one-dimensional Coulomb gas at inverse temperature 
/3. iDumitriu and EdelmanI ( 12002| ) discovered a that ([T]) is the eigenvalue distribution for the 
tridiagonal matrix 
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where the random variables gi,g2, - ■ ■ ,gn are independent Gaussians with mean and vari- 
ance 2 and X/?, X213, ■ ■ ■ , X(n-i)i3 are independents x random variables indexed by the shape 
parameter. 

When n t 00, the largest eigenvalue centered by and sc aled by n^^^ conve r ges in 
law to the Tracy-Widoni(/3) d istribution. This was first shown in iTracy and WidomI (jl994j ) 
a ndlTracv and WidomI (19961) for the cases /3 = 1, 2 or 4, where exact formulae are available. 



Ramirez. Rider and Virag] (120061 ) extended this result for all the (3. They show that the 
rescaled rescaled operator: 



n 



1/6 



converges to the stochastic Airy operator {SAEjs 



(2) 



in the appropriate sense (here B' is a white noise). In particular, the low-lying eigenvalues 
of converge in law to those of Hp. Thanks to the Ricatti transform, the eigenvalues 
of SAEp can be reinter preted in terrn s of th e explosion probabilities of a one-dimensional 



diffusion. In particular. 



Ramirez et al. 



jaooel) 



show that 



f>{TWp>a) = P+oo(X blows up in a finite time). 



(3) 



where X is the diffusion 



dX{t) = {t + a- X\t))dt + -^dBit), 
X(0) = 00. 



(4) 



Note also that X — -^B satisfies an ODE, simulated on the front page with (3 = 00 and 
2. The starting time of the separatrix is distributed as — TW^. 
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)• 



Valko and Virad fl201ol ). Ipvsonl fll962h and 



In this article, we study the diffusion (j4]) in order to obtain the right tail of the Tracy- 
Widom law. Our main tool will be the Cameron-Martin-Girsanov theorem: it permits us to 
change the drift coefficient of the diffusion and evaluate the probability of explosion using 
the new process. 

Using th e variational characte rization of the eigenvalues of SAE^ ([2]) and an analysis of 
the SDK 



Ramirez et al. 



( I2OO6I ) show that as a — )• 00 we have 



P {TWp < -a) 
P {TWp > a) 



exp (-^/3a^(l + 0(1))) , and 
exp(-|M2(l + o(l))). 



(5) 



While we were finishing this article, iBorot. Eynard. Majumdar and Nadall (l2010[ l , in a physics 
paper, using completely different methods, calculated more precise asymptotics for the left 
tail of the Tracy- Widom distribution. 

In this paper we evaluate the exponent of the polynomial factor in the asymptotics of 
the right tail. 



Theorem 1. When a — t- +00, we have 



P {TWfs >a) = a-3^/^exp (-^/^a'^' + 0(v^)^ 



(6) 



This generalizes, in a less precise form, a result that follows fro m Painleve asymptotics 
for the case (3 = 2 (see the slide 3 of the presentation of iBaikI (120061 )). 

PiTW, >a) = ^ exp (^~a'/' + O {a-'/')^ . 
The structure of the proof of Theorem [1] is contained in Section [2l 



Preliminaries and notation. For every initial condition in [—00, +00], the SDK (jlj) ad- 
mits a unique s olutio n, and this solution is increasing in a for each time t (see Fact 3.1 in 



Ramirez et al. 



(I2OO6I )). From now on, we denote by (fi, J^, P(t .j.)) the probability space on 
which the solution of this SDK X begins at time t with the value Xt = x almost surely, and 
^[t,x) its corresponding expectation {x G [—00, +00]). When the starting time is t = 0, we 
simply write and E^;. 

The first passage time to a level x G [—00,00] for the diffusion X will be denoted : = 
inf{s > 0, X, = x}. 
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Throughout this paper, we study many solutions of stochastic differential equations by 
comparing them to expressions involving Brownian motion. The letter B will denote a 
standard Brownian motion on the probability space {Q,J^,P). We will use the following 
easy estimates. For the upper bounds, these inequalities hold for every x > 0: 



< e 2 



P{Bi > x) 

P (suPte[o,i] \Bt\ >x) < 4e-5^'. 
For the lower bounds, there exists Cfem > such that for every e G (0, 1): 



(7) 



P I sup \Bt\ < s \ > exp ( — Cfr 

te[o,i] 



bm n 



In the sequel, asymptotic notation always refers to a — )■ oo unless stated otherwise. Inequal- 
ities are meant to hold for all large enough a. 

2 Proof of the Theorem [1] 



parabola 
where 
the drift 
cancels 




Stretch III 



Figure 2: The critical parabola {(t, x) : t + a — = 0} 

This section gives the structure of proof of the main theorem. The proof of technical 
points will be treated in the following sections in chronological order. 

We rely on the characterization Q, and separate our study of the diffusion (jlj) into three 
distinguished parts demarcated by the critical parabola 



{(t,x) -.t + a-x^ = Q} 
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where the drift vanishes (see Figure [2]). The exponential leading term of the asymptotic 
comes from the part inside the parabola (Stretch II): the drift is positive and makes it 
difficult for the particle to go down. One part of the logarithmic term comes from the time 
it takes to reach the upper part of the parabola (Stretch I): the t-term of the drift adds this 
cost. 

2.1 Upper bound, above the parabola 

At first, let us approximate the critical parabola by the two horizontal lines ^/a and — a/o (as 
the blow-down times will be typically very small). Moreover, the part below the parabola 
gives no contribution for the upper bound, and we use 

Foo(T_oo < +00) < Poo(T_^ < +00). 

The first step is to control the time it takes to reach i/a. Indeed, as the cost for crossing the 
interval [—^/a, ^/a\ increases with time, we need to find a good lower bound for this time. 
A comparison with the solution of an ODE linked to our SDE enables us to have a quite 
precise information: its typical value is 3/8lna/ ^/a, which does not depend on the factor /3. 
It is very unlikely to happen in time faster than 

r_ = (3/8 - 1/ Vh^) In a/ y/a. 



the most probable path of the coiiditionned distribution 




3/81n(a)/v^ «ln(a)/^ 
Figure 3: The most probable path of the conditioned diffusion for a = 100 
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This is the content of Proposition [2l Therefore, using this proposition, the decreasing prop- 
erty in t and the Markov property, we can write: 

Poo {T^ < r-, < oo) < exp (^-l/Je^v^^ ^M^-V^ < (9) 

The asymptotic formula ( !T5|) given by Lemma [6] will highlight the fact that even if the 
process is considered to start immediately at y/a in line (|9]), the award is small (of the order 
exp(0(lna)) compared to the cost it takes to go down quickly. Consequently, with a much 
more significant probability, it will take a longer time than the one considered in (Q to reach 
y/a. Let us find an upper bound for this case: 

Poo (T_^ < oo, > r_) < P,_,^(T_^ < oo) 

Thanks to the Markov property, the process X under the probability measure Pr_,y^ is 
identically distributed with X defined with the same SDE (jlj) where the variable a is replaced 
by a := a + r_ with the initial condition X(0) = y/a. Observe now ^/a < but it does 
not matter as we will be allowed to reduce the interval [—^/a, ^/a\ a bit without affecting the 
relevant terms in our asymptotics. More precisely, the interval we will study for the middle 
is [—y/a + 6, yfa — 5], where b := ^JXrvaja. 

Let denote the first passage times of X. The inequality a/S^ — 5 < \fa gives: 

^ < oo) = P^(r_^ < oo) < P^_,-(T_^+^ < oo). (10) 

2.2 Preliminary upper bound inside the parabola 

Recall (5 := ^In aja we would like to find the asymptotic of: 

< oo) > P^ (T_^ < oo) . (11) 

The key is Girsanov formula. 

Girsanov formula. To evaluate the cost inside the parabola, we use the Cameron-Martin- 
Girsanov formula which allows us to change the drift coefficient of X and evaluate the relevant 
probability by analyzing the new process. The issue is to find a suitable new drift. The best 
would be to have the one corresponding to the conditional distribution of the diffusion X 
under the event it crosses the critical parabola in a finite time: this would lead to an exact 
formula. As we are not able to do that, we use an approximation of the conditional diffusion. 
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In this direction, we introduce a new SDE in which the drift of X is reversed with a correction 
term given by a function y?: 

dYt = {-a + - t + v{Yt)) dt + -^dBt. 

Let denote the first passage time of the process Y to the level x. 

With this diffusion and under some mild assumptions, the Cameron-Martin Girsanov 
formula gives for every non negative measurable function / and every fixed time t > and 
level / G [0,1]: 

^V-a-i if (^"' « < T-VE+i A t)) = E^-i (/ (k, u < r:^+, A t) exp {Gr^^^^AtiY))^ . 

(12) 

More details about this and the application of the Girsanov formula can be found in Section 
14.21 The exact expression of Gt' (Y) contains (3/4 times 

— y/a-\-l 

plus terms involving the function (p. 

Notice that we can already see the correct coefficient in front of the main term. We 
are now confronted with an expectation over the paths of the diffusion Y. To find a good 
estimate of the exponential martingale, we need to control the first passage time to the level 
—^/a + / and check that the diffusion do not go far above ^/a when this time is finite (in 
order to control the last integral in (fT3|) ). We will at first focus on a preliminary bound, for 
which we do not need any result about the first passage time. The price of this approach is 
that it uses a finer control of the paths which go down. 

Control of the paths. To have a good control of the paths, we examine at first a smaller 
interval than [^/a — 6, — a/o + 6]. On this new interval, the diffusion will go down without 
hitting ^/a with a sufficiently large probability. Indeed, we show that for 5 := ^a/Iu a/ ^ 
we have 

f^a-e{T-V-a+e < Oo) = (1 - o(l))P^_, (t_^+, < OO, T_^+, < T^^./s) • (14) 

This is accomplished by two applications of the strong Markov property. From now on, we 
denote T+ = T^_£/2) ^- = ^-^a+e and A = {T+ > T_}, and have: 

^V^-e{T- < oo,^^) < P^H-f (T- < oo) < P^_|(t^_, < cx))p^_,(t_ < oo). 
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Both inequalities use the fact that the hitting probabihty of any level below the starting 
place is decreasing in the starting time of the diffusion X. Rearranging this formula we get 

P^_,(T_ <oo)< — — . 

Lemma m shows that as a — )■ oo the denominator converges to 1. 

Application of the Girsanov formula. We will now study the right hand side of (fT^ 
P^_£(T_ < oo,^) which is the probability of an event under which the absolute value of 
the diffusion is bounded by y/a. 

In order to find an upper bound for the term f|T3l) with I = e, it would be useful to have 
a bound on the time T_. As we do not have any information about that yet, one idea is to 
choose the function ip such that the coefficient appearing in front of the y/aT_ term becomes 
negative. The function ipi of Section U3] works and it gives an upper bound that is sharp up 
to, but not including, the exponent of the polynomial factor. This is the content of Lemma 
ini As a — )• oo we conclude 

P^_,(T_ <oo)< exp + O(lna)^ , e = ^VhT a/i^a. (15) 

In addition. Lemma |6] also shows that with ^ = C3lna/^/a we have 

P^_5(e < T_^+5 < oo) < exp - /31na^ , S = ^In a/a (16) 

i.e. long times have polynomially smaller probability than what we expect for normal times. 

2.3 Final upper bound inside the parabola 

Decomposition according to the time the process spends near ^/a. Let us introduce 
the last passage time to the level ^/a — 6: 



L := sup{t > : Xt = \fa — 5}, 5 = v^ln a/a. 

and use the temporary notation r = c In a/^a. We can use the less precise result and, 
similarly to the part above the parabola, make a change of variable a := a + r. The strong 
Markov property and the monotonicity gives 

(T„^+5 < oo, L > r) < P(, 
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Since ^/a — 5 > y/d, — i we get the upper bound 

P(o,v/a-e) (T_^a+i <oo)< exp (-2/3/3a3/2 - /31na) 

as long as c (depending on (3) is large enough. The last inequality for some constant c > 

follows from the preliminary bound (fT5|) . This again is polynomially smaller than the 
probability we expect for the main event. 

For finer information about the last passage time, one can divide the paths according to 
the value of this last passage time to formalize the idea the process does not earn a lot when 
it stays near ^/a — 6. 

LclnaJ / k k + 1 \ 

The event in the sum implies that X visits ^/a — 6m the time interval but not later. By 
the strong Markov property for the first visit after time kj^fa in that time interval and 
monotonicity, the sum can be bounded above by 

[1 + clnaj P^_5 (l < T_^+5 < oo^ . 

To complete the picture, we find an upper bound of the process between the times and 

1 / ^/a. This will be possible thanks to a comparison with reflected Brownian motion. Indeed, 
the drift is non-positive above the critical parabola, and so up to time I/a/o the process Xt 
started at \/ a + \fa is stochastically dominated by a/c + 1/ ^/a plus reflected Brownian 
motion. This leads to the very rough estimate: 



^ AZTTT^ sup Xt > CsVa < P sup > (c2 - 2)^^ 



< 



exp (-i(c,-2)V/^). 



If C2 is large enough (precisely if C2 > 2/-\/3a/5 + 2), this event becomes negligible 
compared to the probability to cross the whole parabola. Therefore, we can examine the 
studied probability under the event that Xt is bounded from above by C2y/a for t < 1/y/a. 

We denote by C the event under which the above conditions are satisfied 



C := < L < l/y/a, sup Xt < C2\fa, 

[ iG[0,l/v^ 

We have just seen that 

< oo) < (2clna)P^_5 (T_^+5 < oo, C) + O (exp(-2/3/3a3/2 - /? Ina)) . 
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Application of the Girsanov formula. We will apply the Girsanov formula with a 
function ip = such that it compensates exactly the integral 

The suitable function blows up at ^/a and y/a. Therefore we only use it in the interval 
[—\fa + (5, y/a — 5] and set := outside [— -y/a, a/o]. Partially because of those blowups, 
this function creates error terms involving the first passage time to the level + which, 
if finite, by ( !T6|) can be assumed to satisfy T_^^s < ^. Girsanov's formula applied to the 
event {T^_^_s < ^, C} with (ITS]) leads to the fundamental upper bound of Proposition |H1 

{T-^+s < < exp - ^/31na + 0(v^)^ . (19) 

Conclusion for the upper bound. Using f|T9|) with d as in the inequality (fTOj) of the 

part above the parabola, we deduce the upper bound part of IQ. 

2.4 Outline of the lower bound 

The lower bound, as often in the literature, is easier. It suffices to consider the most probable 
paths. For the part above the parabola, we can write the inequality 

P»(T_ < CO) > (t^ < P(_^^) (r_ < CO) 

and use Proposition |2] to bound the first factor. The second factor can be bounded below by 
the following: 

(^v^-^ < X ^V-.-s {t_^,Vs < e) X P(,>VI+^)(^-oo < oo). (20) 

where d := 3/8 In a/ y/a + 1/y/a. 

A domination by a Brownian motion with drift permits to deal with the first term of fl20l) 
which is of the order exp(— 0(Vln a)). The middle term can be controlled with the same 
event C introduced for the upper bound with a replaced by d. We apply Girsanov formula 
directly with the SDE used for the precise result of the upper bound to obtain Proposition 

El 

^V-a^s {T-V-^Vs < ^) > - ^/31na + o(lnS)) P^_,- (^^^^ < |, C') 

(recall the "prime" notation deals with the "new" diffusion Y). 
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A comparison with the solution of a simple differential equation will show that the solution 
of the new SDE indeed has a "large" probability to go down to — + ^ before the time ^. 
This is the content of Lemma [91 

We conclude the proof of the lower bound by checking that the last term of fl20|) is 
also negligible: the proof is similar to the study above the parabola and can be found in 
Proposition (TUi □ 



3 Above the parabola 

We show at first that we need a certain amount of time to reach the level a/o, typically 
a time r := 3/8 In a/ The following Proposition proves indeed that the probability the 
process hits a/o significantly before r is small, but becomes quite large if this hitting happens 
around the time r. 

Proposition 2. The following upper hound holds for all sufficiently large a: 

There exists cq > depending only on /3 such that we also have the lower hound: 

/ 31na\ 1 

Proof. First part. Fix a G (0, 1/8) and let 

3 \ Ina 
r' := I a 



'a 

It helps to remove the time dependence from the drift coefficient of (jl]). In this direction, 
consider the SDE: 

' dY, ={a- Y,')dt + -^dB, ^^^^ 
Yo = +00. 

The process X stochastically dominates Y. Therefore, if is the first passage time to a/o 
for the diffusion Y, we have: 



Poo (T^<r) <P(T^<r'). 



Now, we study the difference Zt := Yt — -^Bt where B is the Brownian motion driving 
F in ([23]). It satisfies the (random) ODE: 



dZ, 



dt. (24) 
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Define M := sup{|i?(|,t G [0,r']} and notice that thanks to the Brownian tail ([7]), 

P(M > 1) < 4exp ^ ^ 



2(3/8 -or) Inay ■ 

By definition, for all t G [0, r' AT^], the process Zt is above i/a — 2/a/5M and conse 
above the (random) solution of the differential equation: 

F'{t) = a- CP it) 
F(0) = +00 

where C has the following expression: 



This differential equation admits almost surely the unique solution 



Wt > 0, F{t) = W^coth(v^t). 



Hence, 



<p(FM-v/^<-^,mf,B,). 

Let us compute F{t'): 
Under the event {M < 1}, 



This implies: 



and 



exp (^-2r' = exp (^-2(3/8 - cr) In a + O (^-^^ ^ ^ 



1 ^^/Ina 
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Taylor expansions ( 126|) and (127|) give: 

Inequality (125|) becomes: 

Poo (T^<r') < P (^F(r') - < -^M, M<1^ + P(M > 1) 



<p(v^-^ + 0(i=) <M, M<l)+P(M>1) 



< 4exp — 13- h o 



3 In a V In a 



If we take a = we have: 

Vina 



and so inequality ( 1211) holds. 

Second part. The proof is similar. Let us check the main lines. At first, we have: 

p» (r^ < > P» (r^,^ ^ ^ - 7h) " (^^ - 7s) ■ 

Now the process {Xt, t G [0,T^]) starting at a value above -y/a has a non-positive drift and 
is therefore stochastically dominated by Brownian motion. Thus the second factor can be 
bounded below by 



(29) 



For the first factor, instead of the SDE fl23|) we choose: 

dYt = (a + r - Y^)dt + -^dBt 
^Yo = +00 

and study the difference 

Set M' := sup{St, t G [0, r - 1/V«]}. For every t G [0, (r - 1/v/a) A T^], the process is 
below 

a + r 



F(t) := J ^ coth (v/(a + r)Ct) 
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where 



C := 1 



M' 



The Taylor expansion of F{t — under {M' < 1} gives: 



Therefore 



P, 



"\/a 



1 \ 



aV4 ^ 



1 



15 



Since M' — B{t — 1/ y/a) has the same law as the reflected Brownian motion at time r — 1/ y/a 
and 15 — 2e^ > 0, we get the lower bound 



2 

7^ 



a 



15-2e2\ 1 
< ^ 1 > co- 



y \/iiia 

Here cq in represents an adequate constant depending only on j3. 



(30) 
□ 



4 Inside the parabola 

The exponential cost comes from this stretch. This section will be devoted to the proof of 
the following Proposition: 

Proposition 3. Recall 6 v^lna/ we have: 

F^_s{T_^+s < oo) = exp (^—Pa^^' -p\na + 0{^/^^)^ . 

4.1 Control of the path behavior 

Here we show a lemma about the return to —y/a + e. 

Lemma 4. For e :— -^yAna/-^ as a ^ oo we have P^_£/2(^i/a-e — oo) 1- 

Proof. Certainly, the probability that X begun at ^/a — e/2 never reaches ^Ja — e is bounded 
below by the same probability where X is replaced by its reflected (downward) at ^Ja — £/2 
version. Further, when restricted to the space interval \^Ja — e, ^Ja — e/2]., the X-diffusion 
has drift everywhere bounded below by t + 1/2^ e. Thus, we may consider instead the 
same probability for the appropriate reflected Brownian motion with quadratic drift. 
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To formalize this, it is convenient to shift orientation. Let now 

X := -^Bit) -\t^-qt. q ■■= Iv^e. 

Let X* denote reflected (upward) at the origin. Namely, 

X*(t)=X(t)-infX(s). (31) 

s<t 

If we can show that P{X* never reaches e/2) tends to 1 when a — )■ oo, then it will also be 
the case for the X-probability in question. 

We need to introduce the first hitting time of level y for the new process X: Ty := mi{t > 
: X{t) = y}. For each n G N, define the event: 

Dn = {X{t) hits — (n — l)£/4 for some t between r^ne/i and T_(^n+i)e/4:} ■ 

From the representation ( ISTl) . one can see that {supi>oX*(t) > e/2} implies that some Dn 
must occur. Indeed, for this X must go above its past minimum by at least e/2, so in this 
case it would have to retreat at least one level before establishing reaching a new minimum 
level among multiples of e/4. Define the event 

£= ^X{t) > - V-2gt-l 

The event S'^ is equivalent to the Brownian motion -^B{t) hitting a line of slope —q starting 
at —1. Since |-B(t)| is sublinear, this will not happen for large enough q. So 

P{£^) when a +oo. (32) 

On 8, when the following term under the square root is positive, we have 



T-x > 72^2^2 + X - 1 - 2q 

Assume that g > 1. A calculation shows that if g < \/x/2 then the above inequality imphes 
'^-x > ^/x/2. So on S for all a; > we have 

+ g > g V a/x/2, 

and so for alH > 

X(r_,. + t)- X{t_,) < -^B{t.x + t)-{qy ^/2)t. (33) 
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Setting X = en/4 we see that for all n > 

P (^Dn n £ I J^T_„^/4^ < P (the process on the right of (!33|) hits e/4 before — e/4) . 

The distribution of that process is just Brownian motion with drift. By a stopping time 
argument for the exponential martingale exp{'yBt — ^7^/2) with 'J = {q ^ ^/nej 4^^f^ the 
above probability equals 1/(1 + e^^^l^) < e-^^l^l^. We get 

P{Pn n ^) < exp (-(g V v^/4)/3£/8) . 

Recall that 



P (sup X\t) > e/2 ) < PiS") + y P(^ n D 



The sum of terms where ne < (4g)^ is bounded above by (1 + (4g)^£ ^)e ^'J^/^. The sum of 
the rest is not more than 



exp (-e3/2/3v^/32) < Ci^ exp (-eg/3/8) 

n>(4(j)2/£ 



where ci depends on /3 only. Replacing e and g by their expressions in terms of a, and using 
(l32l) we obtain: 

P f SUpX*(t) > e/2] < P{S') + a3/4+o(l)e-(/3/16)(16//3)lna _ q 

4.2 Application of the Girsanov formula 

For every ip G C^(M, M) such that sup^gjg |v9(a;)| < y/a (the function ip is in fact small 
compared to the other terms, and it will be chosen after), we would like to consider the 
following SDE (defined on the same probability spaces as (jl])): 

dYt = {-a + - t + ifiYt)) dt + ^dBt (34) 

VP 

Remark 5. The drift of this SDE is the reversal of the drift in the initial SDE The 
solution of the new SDE starting around y/a is a good candidate for the process X con- 
ditioned to blow up to —00 when a goes to +00. Its expression comes from minimizing 
(approximately) the potential: 



{g'{u) — {u + a — g {u)) )du. 
over the set of functions g such that g{Q) = ^/a and g{s) = —^fa. 
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If we look at events under which the diffusion is bounded, it is easy to modify the drift 
of the new diffusion outside the studied domain and prove that the Novikov condition is 
satisfied as long as the examined events are in the space Tt for a fixed t > 0. Let us fix a 
time t > 0, a level / G (0, 1) and denote by T-j. := Tj_(^_;) the first passage times to ±(1/0 — /) 
for the diffusion X (respectively := fo^' We take an event E G J^t^J^T'_ under 

which the paths of the diffusion are bounded by a deterministic value, which can depend on 
a and (3. Since E is J-^- measurable, Girsanov's theorem gives 

The Radon Nikodym density exp{G.MAT'_0^)) is a bounded martingale, and E is J^t'_- 
measurable, so by the optional stopping theorem the above quantity equals 

E^_i(E^_i(lEexpiGtiY))\TT'_)) = ^v^-i exp (Gt'_m{Y))) . (35) 

In the following, we consider E of the form E = {T_ < 00} (1 Ei. The assumption on E 
requires Ei being an event under which the diffusion is bounded from above. Taking the 
limit t — )■ +00 leads to the fundamental formula: 

P^_,(T_ < oo,Ei) = E^_, (ir <oo,i?( exp (Gr (F))) . (36) 

Thanks to Ito's formula we can write the exponential martingale ^Gt^ {¥) as 

2 / {t + a-Y^^)dYt (37) 
Jo 

+ (f){Yo) - (j){Yr) 

rpl 

+ ^ ' -^if'iY,) + ]^ip{Ytf + ipiY,){Y^ -a-t) dt, (38) 

where 99' denotes the derivative of and the indefinite integral. Again by Ito's formula, 
we can compute the first term fl371) above. 



2 /8 \ f - 

(EZD = 2a{Yr - lo) - 3 (>^r " ^'o') + U " j Ytdt + 2rYr. 

Replacing Yt'_ by its value, we obtain the expression (valid under E): 

m = -t^' + + 2ir- 2^T'_ + (I - 2) f^' Ytdt. (39) 
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4.3 The preliminary upper bound 

Let Ci be a constant such that Ci > (|8//3 — 2| — 2)V0. Recall that 



6 := ^Ina/a, e = —=V\na/ \fa 



and take the function defined by 

^r-.x^l-^ (40) 
[0 if a; (-^a, ^a) 

and extend this function on the entire real line such that <y9i remains a smooth function 
supported on [—^0,^0] (this is possible for all large enough "a"). Of course, there are 
many functions satisfying the above conditions but we just need to fix one. Let 0i be an 
antiderivative of 

A first step is to prove a less precise upper bound which does not give us the constant in 
front of the logarithm term of ([6]). In this subsection, the notations F, Ti, will always 
refer to the definitions using this particular 9?!. We have the events 

V w 

Lemma 6. (a) The following inequality holds: 

Pv^-.(T.-v^ < 00, ^) < exp (-^/3a'/' + O(lna)^ . (41) 

(b) for some C3 > 1 we also have 

P^_,(c3lna/v/^ < Ts_^ < 00, C) < exp (-^/Sa'^' - /31na^ . (42) 

Part (a) with Lemma H] immediately give the following. 
Corollary 7. The following upper bound holds: 

P^_.(r.<oo,<exp(-H^«=/=.0(l„.)). 

Proof. Part (a). First let T_ = T^_^. Consider the process Y defined with the function 
ifi. The equahty ( l35l) gives: 

^V-a-e{T- <oo,A)= ^^-,(l{Ti<oo^'} exp (Gri(r)) 

with Gt'A^) given by (EMH]). 



To find an upper bound of tlie last term (1381) in G, we remark tliat tlie cliosen function 
(fi on [—y/a + e, y/a — e/2] attains its maximum at time y/a — e/2. Under the event A', it 
gives: 

'2c2 + 8//3ci 



Moreover: 

Thanks to ([39]), g3]) and (glD, we deduce: 



- Cl^/a T'_ 



ci In ( 2^ + 1 



(43) 



(44) 



^Gri {Y) + ^-a"^ < A^e' + 2eT'_ + 



2 - ci v^t: 



We now take ci such that |8//3 — 2|— 2 — ci<0 and the coefficient in front of y/aT'_ becomes 
negative and dominates the terms involving T'_. The last one creates the logarithmic error, 
and fHTj) follows from 

Gt'_{Y) < -'^I3a"^+ i^ci(3 + lQ\ lna + o(lna). 



Part (b). Now let T_ = T^_^. Just as in part (a), we need to bound the Girsanov terms. 
To find an upper bound of the last term (|38l) . namely 

rpl 

|^'(^*) + ^^(^*)' + ^(^*)(^*' -a-t)dt 

note that ip^ and if' are both uniformly o{y/a). On the other hand, we have Lp{Yt){a — Y"^) 
is nonnegative. Moreover, it is greater than Cl^/a + o(i/a) as long as L' < t < T'_. So on C 
we have the lower bound 



< o(v^)t: -civ^(r: -L'). 



(45) 



The other inequalities are similar to part (a) with 6 replaced by e, except the last term in 
(|39|) gives an extra term due to the fact that Y is only bounded by C2y/a up to time L'. 
Thanks to ([39]), (|45]) and (gl]), we deduce: 



^Gt'_{Y) + ^a'/' <AV^6' + 2ST'_ + 
P "J 



- 2 - ci J v^t: 

+ (ci + C2)V^L' + o{^)r_ + ci In + 1^ 
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From part (a), we have c:= |8//3 — 2|— 2 — Ci < 0. We chose C3 > 1 large enough so that the 
terms involving Ti > C3 In a/a together with the In a term coming from the antiderivative 
are less than — /31na, i.e. f C3C + < —/3. This completes the proof of fj42|) since 

Gt^(Y) < --f3a^/^ - f3\na + o(\na). □ 
3 

4.4 Precise asymptotics for the exponent 

Recall 6 := ^\na/a, L the last passage time to ^/a~6, and the event introduced for technical 
reasons: 

C := < L < l/i/a, sup Xt < C2\fa 
\ te[o,i/v^ 

defined in the outline of the proof. We will study in this section P^_a (T_^^5 < 00, C). 

In order to obtain the coefficient in front of the logarithm term, we need to be more precise 
in our analysis and we will look more carefully at T__, the first passage time to — -\/a + b of 
X. 

Our tool is again the Cameron-Martin-Girsanov formula with a drift containing a different 
function Let us define <y92 in the following way: 

^2:x^\ n H ( r r\ ^^^^ 

y if X ^ (— V a, v«) 

and extend it such that it remains a smooth function on M satisfying: sup |<y9| < ^fa and 
sup Iv^'l < (this is possible for a large enough "a"). Similarly to the previous subsection, 
the notations y, Ti, C etc. refers to definitions with this chosen function. 

Proposition 8. We have 

P^_5 (r_ <C3lna/v^, C) 

/ 2 3 \ 

= exp I - -/3a^/^ - -/31na + O(Vlna) )P^_5 (Ti < C3 In a/^a, C) 

Proof. Let us compute the new Radon-Nikodym derivative according to the position of Y 
using the relations fl37ti38p and (fT3|) . 

At first, the term 0(^0) ~ 4'{Yt'_) is equal to —3/2 In a. Moreover, 



a)^{y) = 0. 
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Consequently, there exists a constant c' > 0, depending only on /3, such that for every 
y e [-y/a + 6,y/a-6], 



< 



c'a 2c' 

2\2 - 



[a-y 



There is another constant c" > such that 



u ip{Yu)du 



< — T'^ 



For every y > ^/a — 5, 



'^^\y) + ]^^{yf + ^{y){y''-a-t) 



< 



/3 



+ 1 Va 



Putting all together and using the upper bound on the last passage time to ^/a — 5 contained 
in C, we obtain: 



If {T'_ < C3 In a/Va} holds, 

2c' 



3 



52 



s 



2 



□ 



— r + -r:^ + 25Ti + -5=^ < 2c' Vh^ + 0(1). 
0"* 3 

Under {T'_ < C3 lna/-\/a} fl C, we conclude 

Gy, (F) = --/3a3/2 - -/31na + f^hT^) 
3 8 V / 

To complete the study inside the parabola for the lower bound, we prove: 
Lemma 9. There exists C4 > depending only on /3 such that with ^ = C3lna/^/a 

^V^-siTL < C) > exp (-C4V^) . 

Proof. For the lower bound we can replace the event C by the event that T| = T'^_^^^ is 
infinite, i.e. the corresponding level is never hit. We will show that this events happens as 
long as 

M:=sup{|Ei|, te[0,C]}<S^/^/5, 

which by the Brownian motion estimate ([H]) has the right probability. 

Let ^ = C3 \na/y/a. Again, we compare our equation to an ODE. The quantity Z = X — B 
on [0, ^] satisfies 

Z' = -{Z - + t-a< Z^ + a+ ^MZ + ^, 
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so let H be the solution of the (random) ODE: 

H'it) = H\t)-C. 



(47) 

if(0) = v^-5, 



where the random constant satisfies 



C = a-^^M + 0{1). 
V p 

By the same argument of comparison as in the Section |3l when M < c6 the diffusion Y is 
under t H{t) + 2/ y/^Bt up to the minimum of ^ and the exit time from [— i/a — 1, a/o]. 
Therefore we will have < ^, Ti < as long as 

H{^) + ^Bf. < -v^ + 5, and sup (h{s) + -^BA < - 5/2. 

VP se[o,c] V yp J 

Since -ff(s) is decreasing in s, the second event is implied by our assumption on M. 
The solution H takes the form: 

, , ftanh f Vet) - fe) /^r _ r 

H{t) = -VCtanh (Vet - arctanh (6) ) = ^ ^ ^, b = ^ ^ . 

^ ^ 6tanh(yCtj-l VC 



When C3 > 1 we have tanh(vot) = 1 + 0[a ^). So we get the asymptotics 

= -v^ + 2M/v^ + o{6), 

and we indeed have 



AM ^4 



H{0 + —B^<-VE+ — < -5 + 0(5). □ 



5 Under the parabola, lower bound 

We will prove: 

Proposition 10. There exists C5 > depending only on (5 such that, 

P_^+5(T_oo < 00) > exp (^-csvlnaj . 

Proof. Using the strong Markov property and the increasing property, we can lower bound 
the left hand side by 



xP 



/ Ina _ r: Vln a \ {T—OO ^ ^X)) . 
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• The first probabifity gives the main cost. Under this event, the process X is stochastically 
dominated by the drifted Brownian motion: 



^ ^2 
VP 



Thus, 



By the reflection principle, this equals 



P (^^^^e >Bi> ^v^^£^ > y^v^ exp (^-^^V]^^ . 

• For the second part, under the studied event, the diffusion (Xt, t > 1/^/a) is stochastically 
dominated by 

1 2 ^ ^ 

+ t + —Bt - Va - £■ 



4v^ 

Thus the studied probability is bounded from below by a constant depending only on (3. 

• For the last part, we need to compare the diffusion with the solution of a simple differential 
equation. Similarly to the previous comparisons, under the event {X{lna/y/a — —^Ja — 
\J\a.al \fa^ TLoo < l\ (3/8 In a/ v^}, the diffusion X is stochastically dominated by 

Git) + 1l\f^Bi where G is the solution of the differential equation: 



and 



Whenever we have 



M= sup l^^l. 

3 In g 1 
'8 V^i 



^gjg 31na, 



2^/;5 ^ 

the function G blows up to —oo at a time smaller than 3/8 In a/ \fa and the diffusion (X^, t € 
[0, 3/8 lna/\/a]) stays under —^fa. Therefore: 

2 v^na' 



2v^ ^ 

which is greater than a constant depending only on /3. It leads to the result. □ 
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